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We study effects of classical magnetic impurities on the Anderson metal-insulator transition nu¬ 
merically. We find that a small concentration of Heisenberg impurities enhances the critical disorder 
amplitude Wc with increasing exchange coupling strength J. The resulting scaling with J is an¬ 
alyzed which supports an anomalous scaling prediction by Wegner due to the combined breaking 
of time-reversal and spin-rotational symmetry. Moreover, we find that the presence of magnetic 
impurities lowers the critical correlation length exponent u and enhances the multifractality param¬ 
eter oq. The new value of u improves the agreement with the value measured in experiments on 
the metal-insulator transition (MIT) in doped semiconductors like phosphor-doped silicon, where 
a finite density of magnetic moments is known to exist in the vicinity of the MIT. The results are 
obtained by a finite-size scaling analysis of the geometric mean of the local density of states which is 
calculated by means of the kernel polynomial method. We establish this combination of numerical 
techniques as a method to obtain critical properties of disordered systems quantitatively. 
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Experimental studies of uncompensated doped semi¬ 
conductors like Sii_a;Pa, (Si:P) show a metal-insulator 
transition (MIT) as function of dopant density x. This is 
one of the most extensively studied cases of a quantum 
phase transition m- The doping increases the carrier 
density and thereby the conductivity, but also creates the 
onsite disorder potential. The random positioning of the 
dopants results in random hopping amplitudes between 
the dopant sites. Moreover, the electron-electron interac¬ 
tion causes spin and charge correlations. Therefore, the 
MIT can neither be described completely by an Anderson 
MIT (AMIT) nig, driven solely by disorder, nor by a 
correlation-driven MIT, the Mott transition [6]. Taking 
into account both correlations and disorder remains an 
open problem of condensed matter theory [71 ig. 

Thermodynamic measurements prove the presence of 
localized magnetic moments in the metallic regime in the 
excess specific heat and a low temperature divergence of 
the magnetic susceptibility [HEIISI [To] . Indications of 
magnetic moments can also be seen in transport mea¬ 
surements, such as the thermoelectrical Seebeck coeffi¬ 
cient mm- These experiments have been interpreted by 
assuming that up to 10 % of all P-atoms contribute local¬ 
ized paramagnetic moments at the MIT [TTTll4| . which 
originate from the localized states in the tails of the im¬ 
purity band. Thus, it is an essential step in understand¬ 
ing the MIT in doped semiconductors to understand the 
influence of magnetic impurtities on the AMIT. 

Since the work of Hikami et al. m it is known that 
weak localization is suppressed in the presence of a fi¬ 
nite concentration of localized magnetic moments, be¬ 
cause the exchange interaction with the spins of the con¬ 
duction electrons breaks time-reversal invariance. The 


breaking of time-reversal symmetry (TRS) is known to 
weaken Anderson localization, so the localization length 
^ is enhanced. In quasi-lD disordered wires this leads 
to an enhancement ^ with /3 = 1 (/3 = 2) when 

TRS is unbroken (broken) and s = 1 (s = 2) when spin- 
rotational symmetry (SRS) is unbroken (broken) [THlfTT] . 
Here, where Mt is the number of transverse 

channels and le is the elastic mean free path [17]. With 
SRS intact, there are independent channels for the local¬ 
ization of up and down spins. Otherwise the spin-up and 
spin-down channels are mixed, and the electrons have 
effectively twice as many channels, which enhances ^ ac¬ 
cordingly. 

Three-dimensional disordered systems are known to 
bear an AMIT. It can be expected that the breaking 
of TRS and SRS shifts the critical disorder Wc towards 
stronger disorder amplitudes W, which is a measure of 
the width of the distribution of the disorder potential. 
Likewise, the critical electron density n-c in doped semi¬ 
conductors is decreased |T8| [19]. The symmetry class 
of the transition is thereby changing from orthogonal to 
unitary m [nj. In the presence of an external mag¬ 
netic field, this change of symmetry class of the conduc¬ 
tion electrons is governed by the parameter X-q = 
where Ib is the magnetic length. Therefore, in analogy, 
the spin scattering rate due to magnetic impurities 
is expected to enter through the symmetry parameter 
Xs = where Lg = is the spin relaxation 

length, De is the diffusion coefficient and ^ is the cor¬ 
relation (localization) length on the metallic (insulating) 
side of the AMIT m- When Xs > 1, the electron spin 
relaxes before it can cover the area limited by ^ and the 
system is in the unitary regime. 
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The crossover at the mobility edge can be studied 
through a scaling ansatz for the conductivity a on the 
metallic side, as done in Ref. [H] for the case of an ex¬ 
ternal magnetic held. Following this approach, using the 
spin scattering rate r“^, yields f {Xs)/{h^). 

As function of the energy difference XE = E — Em to 
the mobility edge Em one then obtains m 

a{T-^) = AE^^-^>f{T-^AE^) . ( 1 ) 


and (fii G [0, 27r]. ai are the Pauli matrices, so the general 
form of the coupling term, Jiai ■ Si, can be written 

Hs = S'^jA cos 9i ^ a\i, a) (i, cr| 

i \ (T—±l 

-|-sin0j ^ exp(zCT</5i) |i, cr) (i, —cr| . (5) 

cr = ±l / 


Simple scaling theory yields (p = 2v [TH], with v be¬ 
ing the critical exponent describing the divergence of ^ 
at the mobility edge (in the 3D orthogonal universality 
class). Wegner argues that while an external magnetic 
held yields pa. = 2^, the additional SRS breaking by 
magnetic impurities rather yields 

V^s = 2zz -|- 3 , (2) 

in a 2nd-order d = 2 -\- £ expansion dH [20] • Thus, a 
numerical analysis of p in the presence of magnetic mo¬ 
ments is called for. The value of pa in a magnetic held has 
been studied in Ref. [U], and they hnd good agreement 
with Lpa = 2j^ within their numerical accuracy for Wc- 
A random magnetic held should yield the same value. 
On the other hand, if a hnite concentration of classical 
magnetic impurities um with spin S is present, the spin 
relaxation rate is hnite, resulting in a shift of the crit¬ 
ical disorder amplitude Wc that can be obtained from 
the scaling ansatz for the conductivity Q. For classical 
spins with a Heisenberg exchange coupling of strength J, 
it follows 

Wc = W° + W°(^^J^ m, ( 3 ) 

where is the critical disorder strength without mag¬ 
netic impurities, De = Vpr/d, vf is the Fermi veloc¬ 
ity, and T~^ is the total elastic scattering rate. Oc is a 
constant representing the smallest length scale, which is 
identical to the lattice spacing here. 

We start from the Anderson model Hamiltonian |23| . 

Ho = t |j,tT) (i,cr| + ^ F* |f,cr) (z,tT| , (4) 

where |i, a) denotes an electron state with spin a located 
at site i of a 3D cubic lattice with N = sites and 
periodic boundary conditions. For the local potential F, 
random values are drawn from a box distribution of width 
W, Vi G [—bF/2, bF/2], while the hopping amplitude t 
between neighboring lattice sites remains constant. 

We add another term to the Hamiltonian, describing 
a local coupling of the conduction electron spin ai to a 
classical spin (two-fluid model) [101HH El] with Sf = 
= 1 and random orientation (Heisenberg-like), given 
by the (polar and azimuth) angles 9i and ipi. The angles 
are drawn uniformly from the intervals cosF € [—1,1] 


We hx the concentration of sites carrying a magnetic mo¬ 
ment to riM = 5%. Note that this is a realistic value for 
real materials like Si:P [IIHIll- Ji is drawn from a binary 
probability distribution, Ji G {AO}, taking a nonzero 
value with probability riM, for which it conforms to the 
exchange coupling strength J. Eq. @ by itself leads to 
a nonmagnetic scattering rate = 2Tr p{£f)/{ 6h). 
The scattering from the magnetic impurities enhances 
the total scattering rate to = Tq"^ -|- 

The spin-resolved local density of states (LDOS) is 

given by Pi,^{E) = YAZi \{iw\kZ - Ek), where \k) 
denotes the eigenstate with eigenenergy of the Hamil¬ 
tonian E[ = Hq + Elg. We use the kernel polynomial 
method (KPM) [25H27j in combination with the Jackson 
kernel which is known to yield optimal results for the 
calculation of the LDOS, as it smoothens Gibb’s oscilla¬ 
tions most efficiently [25]. The KPM expands the target 
function in a series of Chebychev polynomials that are 
only defined on the interval [—1,1], so a rescaling of the 
original spectrum of El is necessary, which we achieve 
by applying a factor a = 24t to all energies, E = aE. 
The Jackson kernel comes with an energy broadening of 
f] = tt/M at the center of the considered interval {E = 0), 
which is rescaled by the same factor, p = ap |2Hj- 

We consider two ensemble averages: The arithmeti¬ 
cally averaged local density of states (ALDOS), 


1 ^ 

Pav{E) = —— Pn{E) , (6) 

® n=l 

which corresponds to the average density of states 
(ADOS) in the thermodynamic limit (large number of 
samples iVs) [22], and the geometrically averaged local 
density of states (GLDOS), 


Ptyp(L^) = exp 



Ns 

^iogp„(F;) 


n—1 


( 7 ) 


which is also known as the typical density of states [30j . 
Here, the index n takes into account both site index i 
and spin a of the conduction electrons. Ns denotes the 
total number of considered local densities Pn{E) (in this 
work, Ns = 8000). Although the LDOS is known to be 
spatially correlated, to save some computation time, we 
take into account p = 32 random lattice sites from each 
disorder realization in the geometric mean. 
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FIG. 1. (color online) Demonstration of the scaling ansatz Q 
at half filling {E — 0) for three different values J and np = 2. 
The error bars correspond to 95 % conhdence. 

In contrast to the ALDOS, the GLDOS is sensitive to 
the localization character of quantum states. It is re¬ 
duced by both increasing disorder and increasing system 
size within the whole energy spectrum |26j . In the ther¬ 
modynamic limit {L —> oo), the GLDOS approaches zero 
in energy regions of localized states, and a positive value 
in the case of extended states m- For finite system 
sizes the GLDOS stays positive, even for perfectly local¬ 
ized states. In Ref. [30], localised states are detected by 
defining a threshold value for the GLDOS, which is ad¬ 
justed to previously known values for the critical disorder 
Wq. For a quantitative analysis of the critical parame¬ 
ters (including Wc), it is necessary to perform a finite-size 
scaling (FSS) analysis of the GLDOS. To this end, we use 
the scaling ansatz 

r(LF,L) M (8) 

for the GLDOS at half filling {E = 0), where T{E) = 
Ptyp(L')/Pav(L')j d = 3, Ip = (Wc — W)jWc is the reduced 
disorder, v is the correlation/localization length exponent 
[5|, and ao is a multifractality parameter 0 |311 155] . 
We neglect the disorder dependence of pav(A'). By fix¬ 
ing the ratio G = L‘^/M = 1, the function E(x,y) 
([^ only depends on the first argument, x = ipL^/''. 
We expand the function F(x) = F(x, Pa^^irG) to order 
np in X, using a power series, F = 
our fit model ^ possesses Np = np + A parameters, 
(Wc, ao, V, Fo, Fi,--- , F„p). The scaling § is demon¬ 
strated in Fig. [T] for np = 2. 

Tab. [T] summarizes the fit results. Gonsidering five 
system sizes L G {10, 20,..., 50} and eight disorder 
strengths W/t G {15,15.5,..., 18.5}, each fit has access 
to A^d = 40 data points [53]. To assess the quality of 
the fit, the statistic and the goodness of fit (GOF) 
probability Q are computed [28]. For every J, we select 
an optimal series expansion order np G {2,3,4} so that 
\Q — 1/2| is minimized, as Q = 0.5 indicates an optimal 

fit [55] . 

Fig-i shows how the fit parameters Wc, ao and n 
depend on the coupling strength J, using different se¬ 
ries expansion orders np. Our result for Wc at J = 
0 (corresponding to the original Anderson model) is 
jt = 16.52 ± 0.17, which agrees with established. 


TABLE I. Fit parameters Wc, ao and v for different exchange 
coupling strengths J with their standard errors, found by ex¬ 
panding the function F(x) to order np. To assess the quality 
of the fit, we provide the yfi statistic and the goodness of fit 
probability Q |28| . 


J/t 

np 

Welt 


ao 



u 


X 

Q 

0.00 

4 

16.52 

± 

0.17 

4.07 

± 

0.04 

1.48 

± 

0.06 

30.5 

0.54 

0.05 

4 

18.12 

± 

0.28 

4.43 

± 

0.08 

1.37 

± 

0.08 

29.0 

0.62 

0.10 

3 

18.19 

± 

0.21 

4.33 

± 

0.06 

1.40 

± 

0.06 

45.5 

0.07 

0.15 

4 

18.56 

± 

0.31 

4.39 

± 

0.09 

1.32 

± 

0.08 

41.4 

0.12 

0.20 

3 

19.02 

± 

0.25 

4.53 

± 

0.07 

1.32 

± 

0.07 

54.1 

0.01 

0.25 

4 

18.92 

± 

0.26 

4.47 

± 

0.08 

1.29 

± 

0.07 

30.6 

0.54 

0.30 

4 

19.42 

± 

0.31 

4.47 

± 

0.06 

1.50 

± 

0.09 

34.6 

0.34 

0.35 

4 

18.47 

± 

0.25 

4.27 

± 

0.07 

1.25 

± 

0.07 

28.0 

0.67 

0.40 

4 

19.39 

± 

0.21 

4.51 

± 

0.06 

1.31 

± 

0.07 

33.5 

0.39 

0.45 

4 

19.15 

± 

0.32 

4.41 

± 

0.09 

1.34 

± 

0.10 

38.1 

0.21 


more precisely measured values [321 [551 [36] like jt = 
16.530(16.524,16.536) [32]. As the coupling strength J 
is increased, the critical disorder Wc approaches a higher 
value of Wc'^/t « 19.4. This tendency is expected, as the 
symmetry class is changing from orthogonal to unitary. 
Another study considering an external magnetic field has 
found a value of /t « 18.35 [2T]. We conclude that 
the additional SRS breaking of the magnetic impurities 
causes a further increment. This can qualitatively be ex¬ 
pected, since the mixing of the spin-up and spin-down 
channels by the SRS breaking enhances the number of 
available spin channels and thereby weakens the local¬ 
ization [HKn]. 

As can be seen in Fig. [^ the value of ao under¬ 
goes a gradual transition to a larger value by tuning 
up the coupling strength J. Remarkably, this value 
is larger than that of a recent study using multifractal 
analysis of the 3D Anderson model in a magnetic field, 
a^ = 4.094(4.087..4.101) [55]. This suggests that the 
additional spin symmetry breaking enhances ao beyond 
the unitary value as obtained when only TRS is broken. 
For J = 0, our value is in agreement with other studies, 
for example a® = 4.048(4.045,4.050) [55] . 

Our result for the localization length exponent v in the 
orthogonal regime (J = 0) agrees within the achieved 
accuracy with established values [55] [51 IMl 137] like 
i/® = 1.571(1.563,1.579) [57]. However, the error ob¬ 
tained within our method is about 8 times larger than 
that obtained within the well-established transfer ma¬ 
trix method m- During preliminary calculations we 
observed that v is increasing when lowering the ratio 
G = L^/M [55]. This is expected since G is propor¬ 
tional to the ratio between the KPM broadening and the 
average level spacing. For G = 0.1 the broadening should 
be of the order of the mean level spacing, while our calcu¬ 
lations for G = 1 could still mix critical and non-critical 
states [55] . Since the computational effort scales inversely 
linear with G, we are forced to make a tradeoff between 
the largest considered system size N = L^, the chosen 
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FIG. 2. (color online) Dependence of the fit parameters Wc 
(a) ao (b) and ly (c) on the exchange coupling J, using differ¬ 


ent series expansion orders np. The dashed horizontals mark 
established values for the pure Anderson model (A, realised 
by our model for J = 0) [32], a model considering an exter¬ 
nal magnetic field (M) , the 3D orthogonal (O) |33[ I37| 

and the 3D unitary (U) universality class [33] . and the exper¬ 
imental value (exp) [2]. For W^^J), the data with minimal 
\Q — iy2| (a) is fitted to 0 using u = 1.571 (33 (see also 
Tab, [njl. The error bars correspond to 95% confidence. 


value of G and the resulting computation time, so we 
have decided to choose G = 1 for this analysis. Note 
that experimental investigations have always yielded val¬ 
ues of v considerably smaller than theoretical predictions, 
partly because they face a similar problem of low energy 
resolution |2|, just like our numerical method does. 

For disordered systems in a magnetic field (3D unitary 
universality class), values iz smaller than that of the 3D 
orthogonal universality class have been reported 


like = 1.437(1.426,1.448) [33] (marked in Fig. [^). 
Note that within the achieved accuracy, our results for iz 
with magnetic impurities (J > 0) are of similar or smaller 
magnitude, and in good agreement with the experimental 
value jz®’'P « 1.3 [2] of real materials in which magnetic 
impurities are known to exist at the MIT. 

The scaling of 144 with J has been analyzed in Fig. 2a 
Eq. (|^ suggests a scaling 114(J) ^ Hence, we use 


IFc(J) = aJP -b & 


(9) 


for the fit, with /r = 2/(^. The fit results are summarized 
in Tab. |nj The best fit (smallest \Q — 1/2|) is found for 
Wegner’s scaling ([^ with e = 1 m- Also the free fit of 
the parameter fj, shows good agreement with this analytic 
prediction. Our results do clearly not support the rela¬ 
tion ip = 2v [HUSH, which results in GOF probabilities 
that are orders of magnitude away from an acceptable 
range (e.g., Q G [0.1,0.9]). This interpretation remains 


TABLE II. Fit results for 114(J). In the top row, /r is a free 
fit parameter. Otherwise, p = 2/(p is fixed to values (shown in 
bold) according to the given analytic formulas for p |22l I43| . 
using either v = 1.571(1.563,1.579) |37| or our own value for 
J = Q,v = 1.48 ±0.06. 


p =. 



a 


M 


h 



Q 

Free 

fit 

3.40 

± 

0.46 

0.27 ±0.09 

16.52 

± 

0.21 

11.3 

0.13 

to 

3 

3.61 

± 

0.34 

0.33 

16.57 

± 

0.19 

12.0 

0.15 

2iz 


4.52 

± 

0.70 

0.64 

16.89 

± 

0.26 

28.4 

4e-04 

2zz -I- 

3 

3.64 

± 

0.35 

0.34 

16.58 

± 

0.19 

12.2 

0.14 

2zz 


4.62 

± 

0.75 

0.67 

16.93 

± 

0.27 

30.9 

le-04 


intact even when using our own value for the localization 
length exponent v = 1.48 ± 0.06 (for J = 0) instead of 
the value iz = 1.571(1.563,1.579) [^. 

To conclude, we have shown numerically how local 
magnetic moments which break TRS and SRS affect the 
metal-insulator transition in the 3D Anderson model. 
We found that the critical exponent v decreases for in¬ 
creasing coupling strength J and determined its value as 
zz® « 1.3 ± 0.1 for 5% magnetic impurities. Within the 
obtained accuracy, this value agrees with experimental 
results obtained from conductivity scaling at the MIT 
in phosphor-doped silicon [5]. We also find the multi- 
fractality parameter a® « 4.4 ± 0.1 when both TRS and 
SRS are broken, a value larger than the unitary value 
when only TRS is broken [33]. We considered the scal¬ 
ing of the critical disorder amplitude 144(J) and con¬ 
firm an analytical prediction by Wegner m- Thus, the 
present investigation may relate the systematically lower 
values for the critical exponent v found in experiments 
to the presence of a finite density of localised magnetic 
moments. We note that magnetic moments are known to 
form due to the local interaction in localised states mm- 
Recently, it has been shown that the interplay between 
the Kondo screening of magnetic moments and Ander¬ 
son localization may result in a novel quantum phase 
transition [22l [SS] [39] . It remains to combine this effect 
with long-range Coulomb interaction in disordered elec¬ 
tron systems [iniiii] in order to achieve a complete un¬ 
derstanding of the experimental results [5]. Further, we 
established a method to obtain critical properties of dis¬ 
ordered electron systems by a finite-size scaling ansatz for 
the geometric mean of the local density of states, which 
enables us to make use of the kernel polynomial method 
to efficiently calculate the LDOS. This method should be 
further developed in order to reach an accuracy compara¬ 
ble to more established methods like the transfer matrix 
method [42] . 

We gratefully acknowledge useful discussions with 
Georges Bouzerar, Ki-Seok Kim, Hyun-Yong Lee and 
Eduardo Mucciolo. This research was supported by the 
World Class University (WCU) program through the Na¬ 
tional Research Foundation of Korea funded by the Min¬ 
istry of Education, Science and Technology (R31-2008- 






























5 


000-10059-0), Division of Advanced Materials Science. 
The numerical calculations have been performed using 
computational resources of the Computational Laborato¬ 
ries for Analysis, Modeling and Visualization (CLAMV), 
Jacobs University Bremen, Germany. 


* d.junK@lacobs-university.de 

^ s .kettemann@j acobs-university. de 

* slevin@phys.sci.osaka-u.ac.jp 

[1] H. von Lohneysen, Adv. Solid State Phys. 40, 143 (2000). 

[2] H. von Lohneysen, Philos. Transactions Royal Soc. A: 
Math. Phys. Eng. Sci. 356, 139 (1998). 

[3] N. Mott, Le J. de Physique Colloques. 37, 301 (1976). 

[41 B. Kramer, and A. MacKinnon, Reports on Prog. Phys. 

56, 1469 (1993). 

[5] F. Evers, and A. Mirlin, Rev. Mod. Phys. 80, 63 (2008). 

[6] N. Mott, Rev. Mod. Phys. 40, 677 (1968). 

[7] P. A. Lee, and T. V. Ramakrishnan, Rev. Mod. Phys. 

57, 287 (1985). 

[8] D. Belitz, and T. Kirkpatrick, Rev. Mod. Phys. 66 
(1994). 

[9] H. Schlager, and H. Lohneysen, Europhys. Lett. 661 
(1997). 

[10] K. Andres, R. N. Bhatt, P. Goalwin, T. M. Rice, and 
R. E. Walstedt, Phys. Rev. B 24, 244 (1981). 

[11] S. Sachdev, Phys. Rev. B 39, 5297 (1989). 

[12] R. N. Bhatt, and T. M. Rice, Phys. Rev. B 23, 1920 
(1981). 

[13] R. N. Bhatt, and P. A. Lee, Phys. Rev. Lett. 48, 344 
(1982). 

[14] M. Milovanovic, S. Sachdev, and R. N. Bhatt, Phys. Rev. 
Lett. 63, 82 (1989). 

[15] S. Hikami, A. Larkin, and Y. Nagaoka, Prog. Theor. 
Phys. 63, 707 (1980). 

[16] S. Kettemann, Phys. Rev. B 62(20), R13282(R) (2000). 

[17] K. Efetov, Adv. Phys. 32, 53 (1983). 

[18] D. Khmel’nitskii, and A. Larkin, Solid State Commun. 
39, 1069 (1981). 

[19] F. J. Wegner, Nucl. Phys. B 270, 1 (1986). 

[20] There exist also 4th-order d = 2 -|- e expansion results |43| 
for yja and ips that constitute large corrections to the 2nd- 
order result for e = 1. 2nd-order m, 4th-order |43| and 
5th-order [44] results have so far shown poor agreement 
with numerical values [5^ 1551157] for the localization 
length exponent iz. However, a new Borel-Pade analy¬ 
sis has recently been suggested |45j that greatly improves 
agreement between the 5th-order 2-l-e expansion [44] and 
numerical values for v. Thus, it should be checked if this 


new Borel-Pade analysis could lead to refined predictions 
of ipa and 1^3 as well. 

[21] T. Drose, M. Batsch, 1. K. Zharekeshev, and B. Kramer, 
Phys. Rev. B 57, 37 (1998). 

[22] S. Kettemann, E. R. Mucciolo, 1. Varga, and K. Slevin, 
Phys. Rev. B 85, 115112 (2012). 

[23] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[24] M. A. Paalanen, J. E. Graebner, R. N. Bhatt, and 
S. Sachdev, Phys. Rev. Lett. 61, 597 (1988). 

[25] A. Weifie, G. Wellein, A. Alvermann, and H. Fehske, Rev. 
Mod. Phys. 78, 275 (2006). 

[26] D. Jung, G. Czycholl, and S. Kettemann, Int. J. Mod. 
Phys. Conf. Ser. 11, 108 (2012). 

[27] D. Jung, and S. Kettemann, AIP Gonf. Proc. 1610, 77 
(2014). 

[28] See supplemental material at [URL will be inserted by 
publisher] for further details. 

[29] In fact, we are using an alternative, more efficient KPM 
algorithm to calculate the ADOS directly |25| . 

[30] G. Schubert, and H. Fehske, in Quantum and Semi- 
classical Percolation and Breakdown in Disordered Solids, 
Lecture Notes in Physics Vol. 762 (Springer, Berlin, Hei¬ 
delberg, 2009). 

[31] Y. Asada, K. Slevin, and T. Ohtsuki, unpublished (2006). 

[32] A. Rodriguez, L. J. Vasquez, K. Slevin, and R. A. Romer, 
Phys. Rev. B 84, 134209 (2011). 

[33] L. Ujfalusi, and 1. Varga, Phys. Rev. B 91, 184206 (2015). 

[34] Depending on J, W and L, La < L may not always be 
fulfilled in our calculations [28] . 

[35] The optimal value of is the number of degrees of free¬ 
dom k = Ad — Ap, where Ap is the number of fit pa¬ 
rameters. 

[36] K. Slevin, and T. Ohtsuki, Phys. Rev. Lett. 78, 4083 
(1997). 

[37] K. Slevin, and T. Ohtsuki, New J. Phys. 16, 015012 
(2014). 

[38] S. Kettemann, E. R. Mucciolo, and 1. Varga, Phys. Rev. 
Lett. 103, 126401 (2009). 

[39] S. Kettemann, M. E. Raikh, Phys. Rev. Lett. 90, 146601 
(2003). 

[40] A. M. FinkeLstein, JETP Lett. 37, 517 (1983); 40, 796 
(1984). 

[41] Y. Harashima, K. Slevin, Phys. Rev. B 89, 205108 
(2014). 

[42] K. Slevin, and T. Ohtsuki, Phys. Rev. Lett. 82, 382 
(1999). 

[43] F. Wegner, Nucl. Phys. B 280, 210 (1987). 

[44] S. Hikami, Prog. Theor. Phys. Supp. 107, 213 (1992). 

[45] Y. Ueoka, and K. Slevin, J. Phys. Soc. Jpn. 83, 084711 
(2014). 



